Variation of Longitudinal Plasma Wavelength under Irradiation and Double 
Resonance in Coupled Josephson Junctions 
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The effect of electromagnetic wave irradiation on the phase dynamics of intrinsic Josephson junc- 
tions in high temperature superconductors is investigated. We predict three novel effects by varia- 
tion of the radiation amplitude and frequency: changing of the longitudinal plasma wavelength at 
parametric resonance; double resonance of the Josephson oscillations with radiation and longitudi- 
nal plasma wave; charging of superconducting layers in the current interval corresponding to the 
Shapiro step. The "bump" structure in IVC recently observed experimentally is demonstrated. We 
also observe ragged Shapiro steps at double resonance. 



Electrical and magnetic properties of intrinsic Joseph- 
son junctions (UJ) in high temperature superconductors 
(HTSC) are strongly nonlinear and determined by their 
phase dynamics. The phase dynamics of IJJ is used to 
explain the mechanism of coherent terahertz radiation 
which is investigated intensively today. [lJ-Q] One of the 
most spectacular indications of the Josephson effect in 
HTSC is locking of the Josephson oscillations of each 
junction to the frequency of external electromagnetic ra- 
diation. This locking leads to appearance of steps in the 
current voltage characteristics (IVC) at quantized volt- 
ages V n = nhuj/2e, called Shapiro steps (SS). Here lo is 
the frequency of the applied signal, h - the Planck con- 
stant, e - the elementary charge; n - integer number. 0,0] 

Another interesting feature of IJJ is a longitudinal 
plasma wave (LPW) propagating along the c axis. 0, [13] 
It follows from the fact that the thickness of the S-layers 
is comparable to the Debye screening length rjy, and so 
there is no complete screening of electric charge in a sep- 
arate S-layer. The frequency of Josephson oscillations 
loj is determined by the voltage in the junction, and 
at loj = 2ujlpw (wlpw is LPW frequency) the para- 
metric resonance (PR) is realized: the Josephson oscil- 
lations excite the LPW by their periodical actions. IVC 



of IJJ demonstrate the multiple branch structure [1 lMl4 1 
and have a breakpoint (BP) and some breakpoint region 
(BPR) in the outermost branch before transition to the 
inner branch. 

As it is known, the one-dimensional models with cou- 
pling between junctions capture the main features of real 
IJJs, like hysteresis and branching of the IVC, and help 
to understand their physics. An interesting and very im- 
portant fact is that the ID models can also be used to 
describe the properties of a parallel array of Josephson 
junctions, which is often considered as a model for long 
Josephson junctions. In particular, in Refslla. Il6l the 
experimental data have demonstrated a series of reso- 
nances in the IVC of the array and were analyzed using 
the discrete sine-Gordon model and an extension of this 



model which includes a capacitive interaction between 
neighboring Josephson junctions. The parametric insta- 
bilities of an one-dimensional parallel array of N identical 
Josephson junctions were predicted by theoretical anal- 
ysis of the discrete sine-Gordon equation (also known as 
the Frenkel-Kontorova model) and observed experimen- 
tally in Ref|17|. In particular, the novel resonant steps re- 
lated to the parametric instability were found in the IVC 
of discrete Josephson ring. It was verified experimentally 
that such steps occur even if there were no vortices in the 
ring. 

In this letter we investigate the effect of electromag- 
netic irradiation on the phase dynamics of IJJs and tem- 
poral oscillations of the electric charge in superconduct- 
ing layers. The results carry with them several important 
messages concerning the variation of the wavelength of 
the LPW, the charging of the S-layers and their effects 
on IVC of the coupled J J. We investigate the influence 
of the external radiation on the PR by: i) increase of 
amplitude of radiation at fixed frequency ; ii) decrease 
of frequency of radiation which brings SS to the para- 
metric resonance region (PRR). An irradiation leads to 
the decrease of the hysteresis in IVC of IJj[3, El, so 
it's expected that the PR point (loj — 2lolpw) would be 
shifted also, and LPW frequency would increased . As 
far as we know such study have not been done yet. 

To investigate the phase dynamics of IJJ we use one- 
dimensional the CCJJ+DC modeljl^] where the gauge- 
invariant phase differences <pi(t) between S-layers I and 
I + 1 in the presence of electromagnetic irradiation is de- 
scribed by the system of equations: 



difi 
dt 



V l -a(V l+1 + V l . 1 -2Vi) 



dVi dtpi 

— = I - sin ipi - j3— + A sin Lot + I no i, 



(1) 



where t is dimensionless time, normalized to the in- 
verse plasma frequency lo~ x , where lo p = y^2eI c /HC, 
(3 = l/\f]T c , /3 c -McCumber parameter, a gives the cou- 
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pling between junctions @, A is the amplitude of the ra- 
diation. To find the IVC of the stack of IJJ we solve 
this system of nonlinear second-order differential equa- 
tions (1) using the fourth order Runge-Kutta method. 
In our simulations we measure the voltage in units of 
Vq = huj p / (2e), the frequency in units of uj p , the bias cur- 
rent I and the amplitude of radiation A in units of I c . To 
study the time dependence of the electric charge in the 
S-layers, we use the Maxwell equation div(eeoE) = Q, 
where e and £rj are relative dielectric and electric con- 
stants. The charge density Qi in the S-layer I is propor- 
tional to the difference between the voltages V\ and V;+i 
in the neighbor insulating layers Qi = Qoa(Vi + i — V/), 
where Qq = £eoVo/r 2 D . We have taken even number of 
junctions in the stack ./V = 10 to escape the additional 
modulations of the electric charge oscillations in the BPR 
which appear in case of odd number of junctions and 
see clear effect of radiation only. Numerical calculations 
have been done for a stack with the coupling parameter 
a = 0.05, dissipation parameter ft — 0.2 and periodic 
boundary conditions. We note that the results are not 
very sensitive to the parameter values and boundary con- 
ditions in their wide region. The details of the model and 
simulation procedure are presented in Refll4 

It is known that in case of single J J the increase in irra- 
diation amplitude A decreases the hysteresis region, i.e., 
it leads to the decrease of the critical current value and 
the increase of the return current Zp.[l8j| For the stack 
of coupled JJ the external radiation leads additionally 
to the series of novel effects related to the parametric 
resonance and creation of the LPW propagating along 
the c axis.0,[T3] We will describe below the changing of 
LPW wavelength, additional resonances around SS and 
double resonance wj = uj = 2ujlpw with increase in the 
amplitude of radiation A. 
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FIG. 1: IVC of a stack with 10 coupled JJ without irradia- 
tion (curve 1) and under radiation with frequency uj — 2 and 
amplitude A — 0.1 (curve 2) and amplitude A — 0.5 (curve 
3). Insets (a) and (b) illustrate the modulation of the charge 
oscillations at A = 0.1 and the "bump" structure in IVC at 
A = 0.5, relatively. 



First we describe the case uj > 2uj lpw, when the SS 
is above the PRR in IVC. Fig. Q] show three IVC of a 
stack with 10 coupled JJ: without irradiation (curve 1) 
and under radiation with uj = 2, A — 0.1 (curve 2) and 
A = 0.5 (curve 3). At u> = PR is characterized by 
breakpoint current Ib p — 0.28 and breakpoint voltage 
Vb p — 11.51, corresponding to the Josephson frequency 
ujj ~ 1.151. [12] With increase in amplitude of irradiation 
the PRR in IVC is shifted up along the voltage axis. 
As we can see, the first Shapiro step is developed on 
the outermost branch of IVC in the hysteresis region at 
V = oj j * N = 20. Dashed line stresses this fact. Arrows 
indicate the position corresponding to the appearance of 
PR in the stack. Hollow arrow indicate the additional 
PR before SS which we will discuss below. 

Inset (a) shows time dependence of the electric charge 
in two superconducting layers (LI and L2) at A = 0.1. 
Instead of monotonic exponential increase of charge ob- 
served in case without radiation, we see exponential but 
modulated growth. The charges on the neighbor layers 
are equal in magnitude and opposite in sign. This is true 
for all adjacent layers and corresponds to the 7r-mode, 
so the 7r-mode is survived under radiation with A = 0.1. 
Results of FFT analysis (not presented here) of the time 
dependence of voltage V(t) in each JJ and charge Q(t) in 
each iS-layer show that this modulation is due to beating 
between the external and Josephson frequencies. 

The irradiation can change the character of the charge- 
time dependence essentially and bring about a "bump" 
structure on the outermost branch of IVC, as shown in 
the inset (b) at oj = 2 and A = 0.5. At these parameters, 
as was mentioned above, the additional PR appears in 
the stack before Shapiro step (shown by hollow arrow in 
Fig. [I]). We call this resonance a radiation related PR 
(rrPR) to distinguish it from a fundamental PR (fPR) 
manifested also without radiation. The bump structure 
in IVC were recently observed experimentally. 0, HH, [22| 
It is likely that the charging of the S-layers appears not 
only at parametric resonance, but also at other types of 
resonances in coupled JJs. [3, 23, 24 [ This likelihood has 
not yet been investigated in detail. 

Fig.[5]demonstrates the effect of the amplitude increas- 
ing at uj — 2 on the wavelength of the fPR LPW. We 
show that before the resonance region in IVC (Fig. [2^) 
the charge in the layers is zero (in noise level). In the 
growing region of the resonance (Fig. [2})) the amplitude 
of the charge oscillations increases exponentially forming 
the LPW with k = tt (A = 2d). At A — 0.14 the wave- 
length of the created LPW at the parametric resonance 
is changing by the external radiation. The charge distri- 
bution along the stack, presented in Fig. [2jc), illustrates 
the wave with A = lOd. At A = 0.23 we fixed additional 
changes of the LPW: A = lOd A = 5d. 

Results of detailed investigations of the irradiation ef- 
fects at uj = 2 in the amplitude range (0, 0.35) are sum- 
marized in Fig. |3] which shows the variation of the LPW 



3 



Before fPR, A=0 



10 



fPR, A=0.15 fPR, A=0.23 



11 < 
9 < 


• 

. 10 

i 




> 8 






7 < 


> 






< 


• 6 


5 < 




< 


. 4 






3 4 


• 


< 


. 2 






|1 < 


• 




-0 

(a) 



01 0.01 

Q, 



01 0.01 -0.0002 0.000: -0.0002 0.0002 

a. (c) Qi (d) Q, 



12 



11.8 



11.6 



11.4 












I 


co=0 



'1.165 
1 1.160 
1.155 

lsT.= 2a> . =1.151 
1.150 J LPW 

I 1.145 



0.28 



_i_ 



_i_ 



_i_ 



0.282 0.284 0.286 
I 



FIG. 2: Demonstration of changing of LPW wavelength with 
increase of the amplitude of radiation. The numbers count 
the layers in the stack. 



wavelength with A. In the fundamental PR we register 
the following transitions of LPW with increase in A: A = 
2d => A = Wd A = 5rf A = 3d => A = 2d. As we will 
discuss below, the increase in A leads to the appearance 
of an additional PR before SS. The increase in A in the 
interval < A < 0.35 also changes the LPW at this ra- 
diation related resonance: A = lOd A = 5d ^ A = 3d. 
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FIG. 4: IVC of a stack with 10 coupled JJ under radiation 
with amplitude A = 0.005 and different frequencies. Thick 
curve (black online) shows IVC without irradiation. Inset 
stresses the coincidence of the curves before Shapiro step. 
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FIG. 3: Demonstration of changing of LPW wavelength 
at fundamental PR (filled squares) and radiation related PR 
with increase of the amplitude of radiation at uj — 2. 

The double resonance condition ujj = lu = 2u>lpw can 
be approached by decreasing of the radiation frequency. 
In Fig. @] we show the IVCs of a stack with 10 coupled 
JJ under radiation with amplitude A = 0.005 and dif- 
ferent frequencies. The numbers near the corresponding 
curves indicate the frequency of external radiation. Thick 
curve (black online) shows IVC without irradiation. Inset 
stresses the coincidence of the all curves before Shapiro 
steps. The SS does not appear at frequency smaller than 
uj = 1.151, because jump to another branch is happened. 

The approaching to the double resonance demonstrates 
an interesting feature of coupled JJ which is absent in 
case of single J J: when external frequency is close enough 
to the PR condition ujj = 2ujlpw, then charge oscilla- 
tions appear in S-layers in current interval correspond- 
ing to the Shapiro step ("charging of SS"). In our case 



FIG. 5: Demonstration of the Shapiro step "charging" : charge 
oscillations in SS current interval. Inset (a) enlarges the 
charge oscillations in the parametric resonance region, insets 
(b) and (c) enlarges in consecutive order the charge oscilla- 
tions in SS region. 



such charging of SS appears starting from to ~ 1.1555, 
while for fPR without radiation ioj — 1.151. The am- 
plitude of oscillations and current interval of charging 
("width of charging") is growing with approaching of 
double resonance. Fig. [5] demonstrates the charging of SS 
at lj = 1.155. Charge oscillations in S'-layers correspond 
to the 7r-mode of LPW. The enlarged parts of charge-time 
dependence are shown in consecutive order in the insets 
(b) and (c) . The inset (a) enlarges the charge oscillations 
in the time interval, correspond to the bias current close 
to the transition to inner branch and demonstrates that 
the fPR is survive at these radiation parameters. It also 
corresponds to the creation of the 7r-mode of LPW. 
Fig. IS] shows the effect of amplitude increasing at 
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FIG. 6: IVC of a stack with 10 JJ at u> = 1.151 and different 
amplitude of radiation. Inset (a) enlarges the part with small 
A. Inset (b) shows the charge-time dependence at A — 0.06. 



uj = 1.151, i.e. at double resonance conditions. Inset 
(a) enlarges the part of figure with small A. In this case 
even small amplitude of radiation leads to the "charg- 
ing" of Shapiro steps. We see it in the inset (b), where 
the charge-time dependence at A — 0.06 together with 
IVC are presented. At double resonance the transition 
to inner branch is directly happened from the SS. In this 
case the Bessel function dependence of the SS width on 
A for coupled JJ is broken and we observe a "ragged" SS. 
To stress this effect, we show in Fig. [7] A-dependence of 
the bias current at the beginning and at the end of SS 
for single JJ (hollow diamond) and stack of 10 JJ (filled 
squares) at u = 1.151 (Fig. [TJa,) and u = 2 (Fig. [7b). 
Double arrows indicate the corresponding SS width. We 
see that at uj = 1.151 the value of SS width is cut off 
("ragged") in compare with the cases of single J J and 
stack at uj = 2, when SS if far from fundamental PR. 




FIG. 7: (a) Demonstration of SS width changing for single JJ 
and stack of IJJ at u = 1.151 (b) the same at uj = 2 

An increase of the radiation amplitude leads to the 
appearance of the PR before SS (rrPR, I > Iss) and 
the creation of LPW with wavelength A = 5d. Charged 
and ragged SS is still conserved. This situation is demon- 
strated in inset (b) of Fig. [6] which shows the charge-time 
dependence together with IVC at A — 0.06. We see two 
charged regions: the rr-PR region and SS region. Ap- 
pearance of rr-PR in the system before SS is reflected in 



IVC by its deformation (appearance of the breakpoint). 
We indicate such deformation by arrows at A = 0.06 and 
A = 0.08. Further increase of radiation's amplitude leads 
to disappearance of the main SS. We note that we observe 
a "charging" of some SS harmonics as well. Detailed de- 
scription of this phenomena and various manifestation of 
the double resonance in coupled JJ will be considered 
somewhere else.(25| 

As summary, we demonstrated a series of novel effects 
of external electromagnetic irradiation on the phase dy- 
namics of coupled Josephson junctions and their current 
voltage characteristics which are absent in case of single 
junction. The results carry with them several important 
messages about the variation of the wavelength of the 
longitudinal plasma wave with changing the amplitude of 
external radiation, the charging of the S-layers in current 
interval corresponding to the Shapiro step and double 
resonance of the Josephson oscillations with irradiation 
and longitudinal plasma wave. These results should be 
taken into account in further research and applications 
concerning behavior of the intrinsic JJ in HTSC under 
external electromagnetic radiation. 

We thank A. E. Botha, P. Seidel and M. Suzuki for 
fruitful discussion of some results of this paper. 
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